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Abstract 

We calculate the life time of a resonance in our recently developed framework 
for a test-particle description of transport processes for states with continuous 
mass spectra. The result differs from the expression commonly used in trans- 
port simulations but agrees with the one derived by Danielewicz and Pratt 
relating the life time to the scattering phase shift. 
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I. INTRODUCTION 



For the understanding of heavy-ion coUisions semi-classical transport theory has become 
an indispensable tool (cf. e.g. [0,0). The theoretical foundation of this approach goes back 
to the pioneering works on non-equilibrium quantum field theory . While former works 
have focused their attention more or less on the quasi-particle regime (see e.g. |]7|-|r3| and 



references therein) the extension of the formalism to off-shell phenomena has become a topic 
of growing interest in the last few years p!^-p2| since it has been realized that the collision 
rates present in high energetic nucleus-nucleus collisions typically are so large that an on- 
shell approximation seems to be inappropriate. In addition, the resonances excited during 
the reaction may have large decay widths. Therefore, a representation of these states by 
stable particles may not be a proper approximation. 

The usual approach to solve a transport equation is the representation of the phase- 
space density by test-particles (concerning off-shell extensions see [^^-|2^). If not only 
asymptotically stable states but also resonances are simulated by test-particles, one has to 
attribute a finite life time to those 'resonance test-particles', i.e. one has to demand that 
they decay after some time of propagation — provided that they had not suffered from a 
collision prior to their decay. Such resonance test-particles usually have arbitrary invariant 
masses (chosen according to their spectral distribution). The question which (in general 
energy- or mass-dependent) life time has to be attributed to the resonance test-particles is 
under present discussion The commonly used recipe is to identify the life time with 

the inverse decay width of the resonance, evaluated for the respective invariant mass. Near 
threshold the width thus becomes small due to the available phase-space. Hence the life 
time — if identified with the inverse width — becomes large. In it was suggested to 
rather calculate the life time from the time delay that the particles suffer which form the 
resonance. This time delay is given by the energy derivative of the phase shift measured for 
the particles which scatter by forming the resonance as an intermediate state. This quantity 
— and therefore also the life time calculated in this way — vanishes near threshold. Hence 
the two expressions for the life time show a completely different behavior as functions of the 
invariant mass of the resonance. 

Recently a novel approach to the test-particle description of transport processes for states 
with a continuous mass spectrum has been presented in [^,^. There the formalism was 
outlined for a system of non-relativistic (asymptotically stable) states subject to elastic 
collisions. In principle, however, there is no obstacle to treat also resonances in the same 
framework and to answer the question about the proper life time for such states. Since 
in principle the creation and decay of particles is a matter of relativistic field theory we 
find it more appropriate to present the formalism here for a relativistic system. One of 
the goals of was the derivation of the (non-relativistic) equations of motion for the 
test-particles. We do not repeat the straightforward generalization for the relativistic case 
here. The resulting relativistic equations of motion are given in pO[ for energy-momentum 
independent self-energies and in for arbitrary self-energies. We note that the way in 
which these equations are obtained in is different from our approach . Nonetheless 



the resulting equations of motion are the same. (This issue is discussed in detail in ||2l[] .) 

In the next Section we generalize the approach of |^ to a relativistic system (of scalar 
bosons). In this way we set the stage for a proper treatment of resonances and automatically 
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also review the most important points of the formahsm outhned in In Sec. |ITT| we 

calculate the life time of resonances for various situations. We summarize our results in 
Sec. 0. An Appendix is added which contains some technical details. 



II. REVIEWING THE DERIVATION OF THE TRANSPORT EQUATION AND 
THE EFFECTIVE PARTICLE NUMBER DENSITY 



We start out from the Kadanoff-Baym equations 0] for relativistic scalar fields 
(cf. e.g. pr3| , |20[] and references therein) 



;-ni-m^)D<(l,l') 



-□1 



dl 
dl 



1) D<(1, 1') + S<(1, 1) 1') 
1) D>(1, 1') + 1) 1') 



where we have introduced the two-point functions without ordering 

iD<{x,y) = {(j){y) (j){x)) , 
iD>{x,y) = {(j){x) (f){y)) 

and the retarded and advanced quantities 

F'''\x,y) = e(xo - yo) [F>{x,y) - F<{x,y)] , 
F'^^ix^y) = eiyo-xo) [F<{x,y) ~ F> {x,y)] 

with F = D,J]. The self-energy is denoted by S. After Wigner transformation 

F{X, p) = JdSi e'P^'FiX + u/2, X - u/2) 



(2.1) 
(2.2) 



(2.3) 
(2.4) 



(2.5) 
(2.6) 



(2.7) 



of all quantities a gradient expansion of the Kadanoff-Baym equations is performed. Next 
one neglects all contributions which are effectively of second or higher order in the derivative 
with respect to the center-of-mass variable X. In this way one ends up with a transport 
equation (cf. ||2^ for details of the relativistic case, the corresponding non-relativistic case 



is discussed in detail in 



m 



ReE'^\ S<] - A [v^r, (p2 - - ReE'-^*)^<] = y^T S< - iJ:<A (2.8) 



where we have introduced the following real-valued quantities: 
two-point functions 

S<^>{X,p)=zD<'>{X,p), 



(2.9) 



spectral function^ (s = 



^Note that the definition of the spectral function sfightly diff^ers from the non-relativistic version 



|21|. In the non-relativistic limit there is the following connection between relativistic and non- 
relativistic quantities: ^rei — ^ Anr/{^rn), Frei — > Fnr- 
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A(X,p) = i (S>(X.p) - S<iX,p)) = (,._ ,„._Re^rg;)V + .P(.Y.rt - P'l"' 
width 

r(x,p) = (s>(x,p) - s<(x,p)) . (2.11) 

A typical example for the width and the corresponding spectral function is depicted in 
Figs. |l] and ^ respectively, for a toy model displaying the main features of the A baryon. 
Width and spectral function start to deviate from zero at threshold (here pion plus nucleon 
mass). The spectral function grows until the nominal resonance mass is reached. The width 
grows beyond that point. The large energy behavior is chosen such that ^/sT finally reaches 
a constant value. (Note that the displayed mass range for F is larger than the one for the 
spectral function.) Details of the calculation for F (and of the corresponding real part of 
the retarded self-energy) are deferred to the Appendix. They form the basis of all the plots 
which exemplify the relevant quantities. 

It is worth noting that the spectral function is normalized 

T—PoA= /^^=l, (2.12) 

-oo 

if the real part of the retarded self-energy is connected to the width (imaginary part) via a 
(in general n-times subtracted) dispersion relation 



^ReS-*(X;po,p) = Trir^/— ^P, - jPr{X;ko,p) (2.13) 

where V denotes the principal value. An example for the real part of the retarded self-energy 
is shown in Fig. ^ (cf. the Appendix for details). 

The generalized Poisson bracket used in ( p. 81) is defined as 



[A, B] = dx,Adp,B - dp, Adx.B - V xA VpB + VpA V xB . (2.14) 
Note that the drift term inherent to any kinetic equation is hidden in 

[p2 -m\B] = -2podx,B - 2p- V xB . (2.15) 



The r.h.s. of (|2.8| ) can be expressed via the common gain and loss terms using the definitions 
( CT and ( Pll) : 



v/^F 5< - zS<^ = 5< - S> . (2.16) 

loss gain 

In the most general case of an interacting multi-particle state away from thermal equilib- 
rium, the transport equation ( |2.8|) is a rather involved integro-differential equation. Within 



4 



the quasi-particle approximation a quite successful method to solve the corresponding trans- 
port equation is to use a test-particle ansatz for the phase-space density . From the drift 
and (generalized) Vlasov terms of the transport equation one deduces the equations of mo- 
tion for the test-particles while the scattering integrals yield the appropriate cross-sections 
for the collisions among the test-particles. For the case at hand, i.e. beyond the quasi-particle 
approximation, it is suggestive to use also a test-particle representation 

~ ^5(3)(f - Si^pl - E^it)) 5^'\p-m) ■ (2.17) 

i 

Note that here in contrast to the quasi-particle approximation the test-particles are allowed 
to have arbitrary energies not connected to their respective three-momentum by any on- 
shell condition. The remaining question is which quantity / should be represented by test- 
particles. 



So far we have reviewed the derivation of the transport equation ( p.8|) as outlined in 



more detail in [PD|. Albeit starting from the same transport equation we differ from the 
approach of [2C,22] in what follows concerning the interpretation and use of this equation. 
In [l^iill the transport equation was merely used as a tool to determine the evolution of 
test-particles between collisions. In the approach presented here and in ||21| we adopt the 
point of view that the transport equation (|2.8|) once solved as exactly as possible provides 
the appropriate propagation, scattering and decays of the test-particles. The differences 



between |^0| and our approach are discussed in more detail in |21| . One difference concerns 
exactly the question which quantity / should be represented by test-particles. 

There is an obvious constraint for such a quantity: the corresponding number of test- 
particles 

Ntp = jd\ J-0^^pof{t,x;p) (2.18) 

has to remain constant during the propagation of the test-particles and also for elastic 
scattering events encoded by S^ast- Only inelastic scatterings and particle decays both 
contained in 

Awast = ^'^'^ ~ ^dast (2.19) 

should change the number of test-particles Ntp. For the following considerations we need an 
important property of the elastic scattering part of the self-energies which is a consequence 
of detailed balance: 



zSl,,5>-zS>,3,5< =0. (2.20) 



(27r)4 

A natural choice for / seems to be 

f{t,x;p) = S<{t,x-p). (2.21) 



Indeed, the full Kadanoff-Baym equations ( |2.1| , |2.2|) conserve ( p. 181) with the choice (p. 21 



as long as only elastic interactions are present. (This was the reason for the choice ( p.21| ) in 
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p0| , |22[] .) However, we no longer deal with the full Kadanoff-Baym equations but with the 



transport equation (p.8|) obtained from the former by gradient expansion. From (|2.8| ) it is 
straightforward to get 



d 

dtJ^'^ J (27r) 

- i /* - /* /(i^ ('^^^' - ^'^) P-^^' 



with 



Even if the self-energies are restricted to their elastic parts, the r.h.s. of ( |2.22| ) does not 
vanish. In contrast, for the choice 

/ = S'< := 2^sV A{l-K) S< (2.24) 

we deduce from (U) (cf. [H) 

Using ( [2.2(JD we find that this time derivative only deviates from zero in the presence of 
inelasticities. Hence we conclude that ( p.24| ) is a proper choice for the test-particle number 
density. Note that the transport equation (|2.8| ) constitutes an effective theory derived from 
the underlying full quantum field theory by gradient expansion. In this approximation 
scheme the parts of ^^(X, p) with fast oscillations in X are neglected. It is not hard to 
imagine that such a neglect necessitates the redefinition of conserved quantities as — in our 
case at hand — the particle number density. In the following we will refer to 5*^ as the 
effective particle number density. 

It is fortunate to rewrite the transport equation in terms of this quantity: 

[p'-m2-ReS-*,^%]--^[v^r,(/_^2_j^ggrct) S< ^ 



5<ff - zS< S> (2.26) 



with 



S> = 2Ab (2.27) 

and the Bose enhancement factor 
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To obtain (|2.27|) we have expressed in terms of A and 5"^ according to ( p.lO|) . 

is expressed in terms of using ( p.24| ). Note that in (global) thermal equilibrium the 

following relations hold: 



> thermal equilibrium 



S<ip) = 2nBipo)Aip) 
S<M = 2nBipo)2yrsT(p)A\p){l-Kip)) 
ii:<{p) = 2nB{po) ^/sT{p) 

Hp) = i + ?^i?(po) 

where we have introduced the Bose distribution ub- 

The transport equation ( p.26|) forms the basis of the following considerations. 



(2.29) 



III. LIFE TIME OF RESONANCES 
A. Vacuum case 



For the most general off-equilibrium situation the transport equation ( |2.26| ) is rather 



complicated. To focus on the aspect relevant for our purpose we restrict ourselves to the 
following scenario: A bunch of resonances (which do not interact with each other) with 
invariant mass ^/s = \fp^ is uniformly distributed in vacuum and their decay time is de- 
termined. This allows for drastic simplifications of the transport equation. First of all, any 
dependence on x vanishes. Second, the gain rate vanishes since no new resonances are 
formed. Finally, the self-energies do not depend on time. Thus we get 

/ 9ReS-\ 1 g(v/ir)p^-m^-ReS-^ _ ^S> 

(3.1) 

Using p.ll| ) (note: zS^ = 0) and ( p.23| ) one ends up with the very simple relation 



~ ^'^^^ = 2^sTAil-K) = Ail-K) ^^-^^ 

from which the life time can be immediately read off: 

T = 2poAil- K). (3.3) 

As an example the life time of the A baryon is given in Fig. |^ (full line) . The expression on 
the r.h.s. can be rewritten by introducing the phase shift 6 via 

tan 5 = - , . (3.4) 

One finds 
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dpo 



(3.5) 



S is displayed in Fig. ^ It shows the typical behavior of a resonant scattering phase shift 
(here of pion-nucleon scattering), i.e. it is small for low invariant masses and rises strongly 
in the vicinity of the resonance peak mass. There it becomes 90°. 

Relation ( p. 5] ) agrees with the one found by Danielewicz and Pratt |^ but is in striking 
disagreement with the one commonly used in simulations of nucleus-nucleus collisions (see 
also the discussion in ||2^), namely 



(3.6) 



(Note that in the rest frame of the resonance the latter expression reduces to 1 /F while the 
r.h.s. of (U) becomes 2y/^A (1 - K).) 

Relation (|3.5| ) can be easily understood if one considers the scattering of a wave packet 
on a potential [^]: At large distances r from the interaction region the scattered wave is 
given by 



[d'k^{k-p}—fiE,,ek)e-'^''' 
J r 



(3.7) 



where ip denotes the wave packet amplitude. For an amplitude sharply peaked near k = p 
one might expand the scattering amplitude / around Ep-^ 



f{Ek) = exp[ln/(Efc)] ^ exp 



\nfiEp) + iE,-Ep)^\nfiE, 



f{Ep) exp 



(Ek — Ep 



d 

^dEr 



^rifiEp) 



(3.8) 



Hence the time delay which the scattered wave packet suffers is given by 



Im 



d 
dEr 



In f{Ep 



d6 



(3.9) 



where we have introduced the phase shift (for a single channel with angular momentum /) 
via 



/ 



1 

P 



sinS Pi(cos9„) 



(3.10) 



Suppose now that two particles scatter by forming a resonance. In this case the time delay 
caused by their scattering can be identified with the life time of the formed resonance. This 
is exactly the essence of relation ( p. 51) . 



^For simplicity we suppress the angular dependence. 
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Both expressions ( |3.3|) and ( p.6|) are displayed in Fig. |. Obviously their shapes are com- 
pletely different from each other. The most striking difference is their respective threshold 
behavior. While r as given in ( |3.3| ) vanishes at threshold, in contrast l/F diverges there as 
the width vanishes. In addition, r is peaked in the vicinity of the nominal resonance mass. 
There r is roughly twice as large as l/F. The latter quantity decreases monotonically in that 
region. In spite of their completely different shapes, it is interesting to realize that there are 
intuitive interpretations for both expressions ( |3.3D and ( ^.61 ). If a resonance with arbitrary 
invariant mass is considered as a "real particle", it is clear that its life time increases near 
threshold due to the limited phase-space available for the decay products. This consider- 
ation suggests formula ( |3.6| ) corresponding to the dashed line in Fig. ^ If, however, the 
resonance is considered as a quantum mechanical transient state the uncertainty principle 
applies. Hence the more the invariant mass of the resonance deviates from its pole mass 
the less time it is allowed to live. Qualitatively this is the essence of ( |3.3|) displayed by the 
full line in Fig. ^. Clearly our derivation given above supports the latter interpretation, 
in spite of the fact that what we are aiming at is the appropriate life time of (resonance) 
test-particles which intuitively are regarded to be closer to real particles than to quantum 
mechanical states. Nonetheless the transport equation ( |2.26| ) derived from the underlying 



quantum field theory clearly demands that the appropriate life time is given by ( |3.3| ) instead 

of (13). 

At this stage some words of clarification are necessary. For a sharp resonance it is well- 
known (e.g. ||2^) that its "life time" is given by 1/Fon-shcii (in its rest frame) which seems to 
be in line with ( p.6|) but in contrast to ( p.3|) which roughly yields twice as much (cf. Fig. ||). 
This reasoning, however, is misleading. In the preceding paragraphs we studied the life time 
r of a resonance with arbitrary invariant mass. In contrast, the well-known "life time" (r) of 
a sharp resonance is a quantity averaged over the allowed energy (or invariant mass) range. 
As we shall see now both expressions ( p.3|) and (|3.6| ), albeit completely different in shape, 
yield the same average value, if the width is sufficiently small. The natural choice for the 
probability (density) to find a resonance with energy po is given by the spectral function 
( p.lOj ). Recalling the normalization (|2.12|) we define 

ir):^J'MlAr. (3.11) 



To evaluate this average for both expressions ( p. 3] ) and ( ^.61 ) we restrict ourselves to the 
case of small width, i.e. a sharply peaked resonance. In this case we can neglect any energy 
dependence of the width and simply evaluate it at the pole mass: 



on— shell 



T{po = Ep,p) (3.12) 



with Ep := i/m^ + fP. For simplicity we also neglect the real part of the retarded self-energy. 
As already pointed out in the following relations hold (translated here to a relativistic 
system) : 

^ TT 6{p^ - m^) sgn(po) for F, ReE"^"' ^ . (3.13) 
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Hence we find for tlie expression ( |3.3| ) 



M=7«2^r^=(l-A-)^^— 5^. (3.14) 



We get tlie same result using 



p(Po)^i P° ^ . (3.15) 

In tlie rest frame of the resonance this reduces to the well-known expression]^ 

(r) = 1/ron-shcll • (3.16) 

Note, however, that this relation provides only an average value of the life time. It does not 
answer the question how long a resonance with arbitrary invariant mass lives. 

Concerning practical applications, i.e. transport simulations with test-particles, the pre- 
vious considerations suggest the following two possibilities for states with a small width: 

1. One might represent S*^ by test-particles having arbitrary invariant masses. For small 
width the likelihood to produce a test-particle with a mass far away from the pole mass 
is of course negligibly small. Therefore basically all test-particles will have similar (but 
not exactly the same) invariant masses. The life time of a test-particle is given by ( p. 31) . 



2. The second possibility is to attribute the on-shell mass to all test-particles. In this 
case the test-particles represent the integrated quantity / c?(pg) S"^ instead of S^^. 
Consequently all test-particles have to decay with the same integrated (averaged) life 
time (r) = 1/ron-sheii. 

Obviously, concerning states with a large width only the first possibility is reasonable. From 
a practical point of view a word of caution is necessary. To realize the first possibility one 
has to make sure that one has chosen enough test-particles to "fill out" the whole spectral 
distribution. Especially, if the test-particles in the tails of the spectral distribution are 
missing one might obtain an average life time which is too high since (cf. Fig. ^ 

Ts^m? ~ S/Fon-shoU ~ 2 (t) . (3.17) 

Hence for states with a small width the second possibility described above seems to be more 
practical. 

It is worth to discuss the role of the factor (1 — K) in the expression for the life time 
(pl). Due to the contribution (cf. (|]2|)) 



^It is interesting to mention that even for the exemplifying case of the A baryon displayed in the 
figures with its quite sizable on-shell width of 120 MeV the average values in (3.14) and ( 3.15D do 
not deviate much from each other: (r) 1.68 fm/c, (l/F) k, 1.55 fm/c, l/Fon-sheii ~ 1.64fm/c. 
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(3.18) 

dp„ ' 

K diverges at threshold, i.e. where F vanishes. This is shown in Fig. ^. We also deduce from 
that figure that (1 — K) remains positive for all values of s. This is a necessary consistency 
condition since the life time and also the effective particle number density ( p. 241 ) should not 
become negative (cf. the corresponding discussion in [^). In Fig. |^ we compare the life time 



r from ( ^.31 ) with the simpler quantity 2^/s A. As can be seen the main effect of (1 — K) on 
r is to change the threshold behavior of the spectral function for the A (cf. ( |Al|JA2|) ) 



A ^ T k^g^ (^y/s ~ {niN + m.^)^ ^ for y/s ^ ttin + m.„ (3.19) 

1 /2 

r ~ ^(1 - i^) ~ (v^ - (m^r + m^)) for v^^m7v + m^. (3.20) 



to 



This explains the kink of r at threshold as compared to the smooth rising of the spectral 
function. In total, however, r and 2^/s A look rather similar, especially if compared to the 
completely different shape of 1/F displayed in Fig. ^. 



B. More than one decay channel 

We now turn to an extension of the formalism applicable when there exists more than one 
decay channel for the resonance. While the totalliie time is given by ( p.3|) it is also important 
to know the partial decay rates (inverse life times) for the various channels. Formally this 
can be accounted for by realizing that =: S^J,^ is given by a sum of the self-energies 
corresponding to the different decay channels: 

zS>, = 5:^S>=5:2v/^r,. (3.21) 

j j 

The latter equality makes sense if there are no production channels, i.e. = 0. Now it is 



straightforward to deduce from ( p.2|) the evolution equation 



5*5< =-E^^e<ff (3.22) 



with the partial life times 



2p^T,,,A{l-K) F 



tot 



F, F. 



Ttot (3.23) 

where we have attributed the index "tot" to all expressions calculated earlier to distinguish 
them from the expressions belonging to a single one of the decay channels. We note in 
passing that the probability for the resonance to decay into channel j (branching ratio) 
agrees with the probability that one would deduce from (|3.6|) : 

1_ 

(3.24) 



''"tot _ _ Ftot 

1 



Tj i tot p- 
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C. In- medium case 



So far, we have restricted our considerations to resonance decays in vacuum neglecting the 
production channels encoded in S*^. If we want to study creation and decay of resonances in a 
medium, we have to take both loss and gain terms, i.e. the full r.h.s. of ( p.26|) into account. 
For simplicity we ignore the in general rather lengthy l.h.s. of ( |2.26| ) which describes the 
propagation of the test-particles according to drift term, Vlasov term, etc. For our purposes 
it is enough to study the simplified evolution equation 

-2^0 9,5.-^= ^^f,;_^^,) S.:^-ffi<g>. (3.25) 

In general, ifl^ and iS^ are connected by detailed balance. For the case of thermal equilib- 
rium this leads to the Kubo-Martin-Schwinger (KMS) boundary condition (cf. e.g. [p!3| ) 



iS<(p) =«S>(p)e-*'o/^. (3.26) 

To simplify the discussion and to get a closer relation to the life time introduced above we 
assume in the following 

iS<(X,p) <iS>(X,p) (3.27) 

(concerning the thermal state we turn from a Bose to a Boltzmann distribution). Using 
( CT ) we find 

zS>(X,p) ^2v/ir(X,p). (3.28) 
The KMS-condition (detailed balance) becomes 

tf:<{p) ^ 2v^r(p) e-P«/^ (3.29) 



while ( |3.25| ) simplifies to 



dtS<=--S<+'^A (3.30) 

T Pq 



where we have replaced b as defined in (|2.28|) by 1. Hence we deduce that the effective loss 
rate is given by 1/r while the gain rate is iE^/po- However, since \^T{p)/po and not 1/r is 
connected to iS^/po by detailed balance ( |3.29| ) one might worry about the correct thermal 
limit for the particle number density. Thermal equilibrium, i.e. a stationary solution for 
is reached if gain and loss terms cancel each other: 

~-S<+'^A^O. (3.31) 
r Po 

Using the definition of the life time ( |3.3|) and the (simplified) KMS-condition ( p. 29 ) we end 
up with 

S< 4v^r (1 - K) e-P»/^ ^ 4yir A" (1 - K) ub (3.32) 
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where we have replaced the Boltzmann by the Bose function in the last step.0 Indeed this is 
the correct thermal limit for 5*^ as already presented in (|2.29|) . Hence we have shown that 
the combination of loss and gain terms on the r.h.s. of the transport equation ( |2.26| ) indeed 
yields the correct thermal limit for the effective particle number density S^. Turning the 
argument around we conclude that attributing the life time (|3.3| ) to a resonance only leads to 
thermodynamically consistent results if one realizes that the appropriate (effective) particle 
number density is given by introduced in ( |2.24D and not by S^. Assuming erroneously 



that the test-particles which decay according to ( |3.3| ) represent one would essentially 
solve an evolution equation like 

dtS< = --S< + '^A (wrong!) (3.33) 
r Po 



instead of (|3.3CI| ). The thermal limit of the (wrong) equation ( p.33|) is of course given by 



S< 4yir A'^ (1 - K) e-P»/^ (wrong!) (3.34) 
in contrast to the correct behavior (cf. (|2.29|) ) 

S< ^2AnB^ 2 A e-""'!^ . (3.35) 



To conclude we have deduced the expression ( p.3| ) for the life time from the transport 
equation ( p.26|) for S*^. In turn we have also shown that starting from the life time ( p.3|) 
for test-particles one has to realize that these test-particles represent S'^ and not 5^.0 



IV. SUMMARY 

After a short review of the derivation of a transport equation from the underlying quan- 
tum field theoretical Kadanoff-Baym equations we have focused on the test-particle repre- 
sentation of the transport process. At present this seems to be the most promising tool to 
solve the transport equation in practice. We have pointed out (repeating the arguments of 
pT| ) that the proper quantity which should be simulated by test-particles, i.e. which counts 
the test-particles per four-momentum and space volume, is given by the effective particle 



number density S"^ introduced in ( p.24|) . Once this identification has been made it was 



straightforward to derive from the transport equation ( 2.26|) the life time of the resonance 



test-particles. For the vacuum case we obtained expression ( p.3| ) which can be rewritten 
as (|3.5|) confirming the result of Danielewicz and Pratt [^]. As can be seen from Fig. ^ 



^Note that an exact treatment of the r.h.s. of ( ^.251) together with the exact KMS-condition ( ^3.261) 
would have led immediately to the last expression of ( |3.32| ). We made the approximations only to 
get closer contact to the expression for the life time. 

^Concerning these detailed balance considerations the use of l/F for the life time would be com- 
patible with . 
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this result for the hfe time strongly differs from the inverse width prescription commonly 
used in transport simulations. If more than one decay channel for the resonance has to be 
considered, the partial life time is given by ( |3.23| ). The latter formula also applies if (in 
a medium) collisional broadening modifies the total width. Note that for the calculation 
of the life times the required spectral function as well as total and partial widths should 
be determined self-consistently for the in-medium case. Finally we have shown that the 
identification of the proper effective particle number density ( ^.241 ) and the relation for the 
life time ( p.3| ) are intimately connected. On the one hand, the life time formula ( |3.3| ) is 
derived from the transport equation ( |2.26 ) for the effective particle number density ( ^.24 ). 



On the other hand, starting out from the formula for the life time and assuming thermal 
equilibrium leads to an expression for the particle number density which agrees with the 
thermodynamic limit of S^. 

In ||2l| it has been shown that the scattering cross sections for off-shell test-particles^ 



have to be modified to account for the fact that these test-particles represent the effective 
particle number density and not S*^. There are rather simple rules for the required 
modification: 

1. For each incoming off-shell test-particle the cross section has to be divided by 

r = 2y/^TA{l-K) (4.1) 

where all quantities refer to the (in-medium) properties of the off-shell test-particle, 
i.e. r is its total (in-medium) width etc. 

2. For each outgoing off-shell test-particle the Pauli blocking/Bose enhancement factor 
is given by (cf. (|]28D) 

In the same way the result for the life time can be understood: In complete analogy to 
the cross section the decay width Tj has to be divided by ( [4.1|) . In that way the inverse of 
that quantity changes from l/Tj to ( |3.23D which is our expression for the life time. Whether 
these new rules for the calculation of cross sections and decay probabilities lead to observable 
effects in nucleus-nucleus collisions remains to be seen. 



^Here "off-shell test-particles" denote test-particles with arbitrary invariant mass as they appear 
in the definition of / in ( p. 17 ). By chance such test-particles even might be on their mass-shell 
— still we would call them "off-shell test-particles". This phrase merely denotes a concept than a 
single particle, namely the concept that "ofF-shell test-particles" fill the generalized phase-space of 
/oitr- momentum and coordinate space (weighted by their spectral distribution). In contrast "on- 
shell test-particles" live in usual phase-space. While "off-shell test-particles" represent / "on-shell 
test-particles" represent / d{pQ)f. These concepts are also discussed at page |lO| (items 1. and 2.). 
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APPENDIX A: A SIMPLE MODEL FOR THE A RESONANCE 

To exemplify the formulae given above we study a toy model for the A resonance where 
spin (and also isospin) is neglected for simplicity. The width is determined by the A ^ N tt 
decay channel: 

r(g) = To fj'^l for ^/s>mN + ^m^ (Al) 



with 



Here the momentum fcrei of vr and N in the rest frame of the decaying A resonance with 
invariant mass \/s is given by 

_ [{s - (mjv + m^ f) {s - (mjv - m^)^)]^^^ 

krel - . [AS) 



In ([A2D we have included a Blatt-Weisskopf function |P3 parametrized by the length pa- 



rameter R. We adopt the following values to simulate the vacuum decay properties of the A 
resonance: tun = 940 MeV, = 140 MeV, tua = 1232 MeV, To = 120 MeV, and i? = 1 fm. 

Since ^/s^ becomes constant for large energies the real part of the retarded self-energy can 
be calculated by a one-time subtracted dispersion relation according to (|2.13| ). We choose 



the subtraction constant such that the real part of the self-energy vanishes for s = m\. In 
our simple vacuum model the energy integration in (|2.13|) can be rewritten in terms of the 
invariant mass squared. One finally gets: 

ReS-(^) = -{s -ml)V J , . (A4) 

J TV [s — s'j(m^ — s') 

{mjv+m,r)^ 
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FIGURES 
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FIG. 2. Spectral function A as given in ( ^.10[) as a function of the invariant mass. 
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FIG. 3. Real part of the self-energy as given in (|A4|) as a function of the invariant mass. 
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FIG. 4. Life time r (full line) as given in ( |3.3| ) and 1 /F (dashed line) as functions of the invariant 
mass. Note that r is evaluated in the rest frame of the decaying resonance, i.e. po = ^/s. 
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FIG. 7. Life time r (full line) as given in ( |3.3| ) and 2^/sA (dashed line) as functions of the 
invariant mass. Note that r is evaluated in the rest frame of the decaying resonance, i.e. po = ^/s. 
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